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I> ' Abstract 

oo : 

. A newly developed isochrone synthesis algorithm for the photometric evolution of 

\ galaxies is described. Two initial mass functions, IMFs, in particular, the recent IMF 

determined by Kroupa, Tout, and Gilmore, three photometric transformations, and a 
1-Gyr-burst star formation rate, SFR, are used to compute the B — V and V — K color 
index evolution. Non-negligible differences are observed among model results. 
■ In the framework of the galaxy count model by Cohn, Schramm, and Peimbert a 

i ■ simple merging scenario is considered to account for the excess of galaxies observed in 

the blue band counts. The excess is explained by the number and luminosity evolution 
of a group of galaxies called interacting, I. It is assumed that the number of I galaxies 
increases as (1 + z) v due to mergers. Moreover, it is proposed that their characteristic 
. luminosity increases as (1 + z) 3 due to starbursts driven by galaxy-galaxy collision and 

decreases as (1 + z)~ v due to the change in the size of the galaxies. Not much number 
evolution is needed to account for the excess; for example, a model with n = 4.0 predicts 
that about 17 % of the galaxies at z = 0.4 are interacting. Number evolution models 
with a rather high value of n fit better the data; in particular, the model with n = 4.0 
predicts that about 13 % of the galaxies have z > 0.7 in the 21.0 < m&j < 22.5 interval, 
this contrasts with the upper bound of 5 % obtained with the sample of 78 galaxies 
by Colless et al. The excess of high redshift galaxies can not be simply explained 
by changing reasonably the parameters of the luminosity function of I galaxies. This 
result could indicate that mergers are not the whole story. Our best-fit model produces 
the following values for the parameters of the local luminosity function of galaxies: 
a = -1.20, M* = -20.7, and </>* = 1.66 x lO^M^c" 3 (h = 0.5). 



1 INTRODUCTION 

The excess number density of galaxies observed in the blue band (e.g. Tyson 1988; Metcalfe 
et al. 1991; Lilly, Cowie, & Gardner 1991) over the expected from a non-evolving model has 
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generated great interest. Several possible explanations have been suggested to account for 
it: from those that invoke a nonzero cosmo logical constant that provides a greater volume 
to accommodate the excess number of galaxies (Yoshii & Takahara 1988; Fukugita et al. 
1990; Yoshii 1993), to those that assume evolution of the luminosity function of galaxies 
(LFG) via an increase in the number density of galaxies due to mergers (Guiderdoni & 
Rocca-Volmerange 1987; Broadhurst, Ellis & Glazebrook 1992; Colin, Schramm, & Peimbert 
1994, hereafter CSP) or via an increase in the characteristic luminosity of the LFG due to 
starbursts driven possibly by galaxy-galaxy collisions (Carlberg & Chariot 1992). Moreover, 
faint redshift surveys by Broadhurst et al. (1988) and Colless et al. (1990, 1993) place 
galaxies quite near, the median redshift is just ~ 0.2 and ~ 0.3, respectively. In addition, 
no galaxies were found with z > 0.7 brighter than bj = 22.5 in the sample by Colless et al. 
This result imposes a strong constraint on evolving models. In these models a non-negligible 
fraction of galaxies with z > 0.7 is expected, indeed larger than the one predicted by a 
non-evolving model. 

With the advent of infrared arrays and infrared CCD detectors it is also possible to 
count galaxies deeply in the K-hand (Cowie et al. 1993; Soifer et al. 1994). Surprisingly, 
a non-evolving model fits relatively well the K-counts. This almost rules out models with 
nonstandard cosmologies, as the ones with a non-zero cosmological constant, since in this 
case they overestimate the number of galaxies. Yet, the observational restrictions admit one 
more explanation, namely a new population of blue galaxies at z > 0.2 that faded or were 
disrupted so as to become invisible at present (Cowie, Songaila, & Hu 1991; Babul & Rees 
1992). 

Local low surface brightness galaxies have similar physical properties to those observed 
in the faint blue excess population (McGaugh 1994). McGaugh suggests that the excess is 
made up of low surface brightness galaxies and his idea rests on a bivariate LFG; that is, a 
LFG that, in addition to type, depends on surface brightness. 

The result that faint blue galaxies (FBGs) are weakly clustered (Efstathiou et al. 1991) 
has been sometimes taken as an evidence against the merger-driven explanation for the 
observed excess of FBGs (see, for instance, Babul & Rees 1992; McGaugh 1994); the low 
value for the angular correlation function at 30" is explained by Efstathiou et al. assuming 
that most FBGs belong to a population that is weakly clustered and intrinsically faint at 
the present epoch. This result is challenged by a recent work by Cole et al. (1994) where 
no evidence is found for the evolution of the comoving correlation length with redshift to 
B = 22, however a more recent paper by Infante & Pritchet (1994) supports the earlier 
results by Efstathiou et al. 

The hypothesis that the observed excess of FBGs is due to an increase in the number 
density of galaxies that are strongly interacting, the merger-driven excess, is retaken here. A 
phenomenological model was developed in CSP to account for three different observational 
restrictions, namely: (1) the S-counts, (2) the K-counts, and (3) a S-redshift distribution. 
It was due to the increase in the number density of a type of galaxies, called interacting (I), 
that the observations were fitted. 

This kind of model, like others of similar nature, is supported by a recent work on 
high-resolution imaging of faint blue galaxies (Colless et al. 1994) and by HST studies of 
morphology of high redshift galaxies (Griffiths et al. 1994a, b). Colless et al. find, from a 
sample of 26 galaxies which belong to the redshift survey of Colless et al. (1993), 17 galaxies 
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that have an enhanced star formation rate indicated by [Oil] equivalent widths greater than 
20A. About 30% of these galaxies have companions at projected distances closer than 10 h _1 
kpc. 

This model differs from the one by CSP in the following: (1) a different algorithm for 
the computation of the 'passive' luminosity evolution correction is applied, and (2) although 
the total luminosity of the I galaxies, Lj, here is also not conserved, it is not due to that we 
do not take into account the change in size of the I galaxies, as it occurs in CSP, but due to 
the starsbursts driven by galaxy-galaxy collisions. This is what we call 'active' luminosity 
evolution. The model in this paper also differs from the one by CSP in the assumed values 
for the parameters of the luminosity functions of the different types of galaxies; in particular, 
those values assumed for the characteristic luminosities of Sa-c and E/SO galaxies. Because 
they are fainter in CSP, they produce a better fit for the redshift distribution than the one 
we show in this paper. 

An outline of the paper follows: in §2 we briefly discuss our revised color evolution code 
with an isochrone synthesis algorithm, in §3 our merger-driven model is described, and in §4 
we present our conclusions. 

2 COLOR EVOLUTION OF GALAXIES 

2.1 Model 

Our photometric evolution of galaxies code is fully discussed in CSP therefore here just a brief 
review will be given. Unlike CSP where a standard procedure was used (e.g. Bruzual 1983) 
here an isochrone synthesis algorithm is implemented. The advantages of this algorithm 
have been discussed elsewhere (e.g. Chariot & Bruzual 1991), we just would like to mention 
the usefulness of the algorithm in computing color evolution properties for short (~ 10 s yr) 
time-scales of the star formation rate. 

2.1.1 Ingredients 

The ingredients necessary in any spectrophotometric evolution of galaxies code are four: (1) 
a library of stellar evolutionary tracks, (2) a library of stellar spectra, (3) an initial mass 
function (IMF), and (4) a star formation rate (SFR). Our compiled tracks are from Schaller 
et al. (1992) for solar metallicity. These tracks are incomplete and therefore we had to 
complete them. Tracks from other authors were used and various useful simplifications were 
made; in particular, the lifetime of the asymptotic giant branch was left as a free parameter. 
This will be remedied as soon as the evolutionary models for the horizontal and asymptotic 
giant branch from the group of Maeder and collaborators are incorporated. As long as we 
are only interested in galaxy color evolution there is no need for having a library of stellar 
spectra, we just need a photometric transformation from the HR theoretical diagram to the 
observational one. The compilation work by Schmidt-Kaler (1982) has been used now to 
pass from one diagram to the other. We still continue to use the work by Johnson (1966) to 
get the rest of the color indexes. 
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The IMF is normalized so that, 

PTTlup 

/ m<j)(m)dm = 1. (1) 

Although a value of 0.1 M Q is often used for rrii ow , as we do, we should keep in mind that 
the contribution by substellar objects to the integrated infrared light of a galaxy might not 
be negligible. Of course, this depends on the luminosity function of these objects of which 
we know nothing. From chemical evolution considerations Peimbert et al. (1994) limit the 
amount of substellar objects in the solar neighborhood to < 0.02 M Q pc~ 3 . The recent IMF 
determination by Kroupa, Tout, and Gilmore (1993), KTG IMF, for the solar vicinity is used, 
as a good guess, for the IMF for a variety of types of galaxies. Nevertheless, the "standard" 
Salpeter IMF, S IMF, with the exponent x — 1.35 is used for comparison purposes. A value 
of 60.0 M is taken for m up . 

In view of the great uncertainties about stellar formation histories of galaxies, an expo- 
nential law for the SFR is usually taken to model the color evolution of the different types 
of galaxies. The maximum virtue of this scheme is that almost all present galaxy colors can 
be reproduced by varying the time scale, r, of the SFR: from an elliptical galaxy, r = 0.5, 
to an irregular one, r = oo (e.g. Chariot & Bruzual 1991; CSP). These galaxy colors are, 
of course, metallicity-dependent. This effect is now being incorporated in the latest models 
for spectral evolution of galaxies; in particular, for elliptical galaxies see Worthey (1994 and 
references therein). The evolutionary correction, e-correction, to the counts should not worry 
us too much because as we know from CSP there are only three types of galaxies that play 
an important role on the counts. According to our count-scheme, the luminosity of a galaxy 
belonging to the interacting group of galaxies does not evolve; i.e., it has no e-correction. 
We can be confident that the other two types, E/SO and Sa-c, are being modeled relatively 
well; the first, with a burst of star formation that lasts 1 Gyr, the second, with a constant 
SFR. 

2.1.2 Isochrone synthesis algorithm 

The luminosity of a single stellar population (SSP) is given by 

lA\(t) = I™™* 10- 0A ^- M ^(f>(m)dm, (2) 

where m min is the lower limit of the IMF and m max is the maximum mass in the isochrone 
with age t. The integrated luminosity for a composite stellar population, with a history of 
star formation given by ift(t), can be obtained with the convolution integral 

W*)= f i>(t-t')l^{t')dt'. (3) 
Jo 

The isochrone at time t from a complete set of m evolutionary tracks is built as follows: let 
us assume that the tracks are divided into n stages, from zero age main sequence (ZAMS) 
to post-asymptotic giant branch (P-AGB). The indexes i and j number the stages and the 
tracks, respectively, with i = 1, n and j — 1, ...,m. The mass of the star at the stage i is 
given by 

fogrmit) = Atj hg(M j+1 ) + (1 - Aij) log(Mj), (4) 



4 



where 

log^-logt 
,J log -log t j+lti 

In equation (5) tj ti denotes the age of the star of ZAMS-mass Mj at the i th evolutionary 
stage, where 

tj+l,i — t < 

and 

Mj < mi(t) < M j+1 . 

Mi and M m denote the lower and upper limit of the IMF, respectively. This procedure 
is performed for each stage (value of i). The HR diagram physical parameters, logL and 
\ogT e jf, associated to the star of mass mj are obtained by interpolating the tracks between 
the stars of masses Mj and M J+1 . 

To obtain the integrated properties of this SSP we proceed as Bruzual (1992): the number 
of stars of mass rrii, n(rrii), is found by integrating the IMF from mT to m + , where mr = 
^/mi-iTUi and m + = y/m^m^. The luminosity, Iax,i, is then computed by assigning N^j = 
Aijn(mi) stars to the HR parameters of the star of mass Mj+i at the i th stage, and N it j + i = 
(1 — A it j)n(mi) to the star of mass Mj at the same i th stage. We finally obtain Zaa(^) by 
summing over all i's. 

2.2 Results 

In doing photometric evolution of galaxies special care should be put on the photometric 
transformation. This is stressed in Figure 1 where the B-V and V-K evolution is plotted for 
a 1-Gyr-burst SFR, for three different calibrations. The dotted line uses a compilation work 
by Schmidt-Kaler (1982) and Johnson (1966), the solid line utilizes one from Schmidt-Kaler 
(1982) and Bessel and Brett (1988), and the dashed line comes from a compiled calibration 
work by Bruzual and Chariot (1993). The greatest difference among them amounts to ~ 0.25 
mag in V-K at around 1.7 x 10 7 yr; this difference is due to the V-K values adopted for the 
giant luminosity class at the low-temperature regime. The peak at 1.7 x 10 7 yr is produced by 
the supergiant red branch contribution. Hereafter the calibration by Schmidt-Kaler (1982) 
and Bessel and Brett (1988) will be used. In Figure 2 we have plotted the B-V and V-K color 
evolution for a 1-Gyr-burst SFR, in this case for two different IMFs: the KTG IMF (solid 
line) and the Salpeter one (dotted line). The V-K color is slightly bluer when computed 
with the KTG IMF, the greatest difference at ~ 10 8 yr amounts to 0.1 mag. A much smaller 
difference is observed for the case of B-V. 



3 A Merger Model for Galaxy Counts 

This paragraph describes briefly the Qq = 1 merger model for galaxy counts of CSP. Galaxies 
were divided into five different types: E/SO, dE/dSO (spheroidal dwarfs), Sa-c (spirals), dl 
(dwarf irregulars), and I (interacting galaxies). A Schechter analytical representation for the 
galaxy luminosity function was used for all types. It was shown that a gaussian representation 
for the luminosity function (LF) of the brighter galaxies provided a poor fit to the redshift 
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Figure 1: Evolution of the color indices B-V and V-K for a lGyr-burst SFR with three 
different photometric calibrations: Bruzal & Chariot (dashed line), Schmidth-Kaler and 
Bessel & Brett (solid line), and Schmidth-Kaler and Johnson (dotted line). 
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Figure 2: As in Fig. 1 we have plotted the evolution of the color indices B-V and V-K for 
a 1-Gyr-burst SFR. The two curves come from two different IMFs: KTG IMF (solid line) 
and Salpeter IMF (dotted line). The calibration from Smichdt-Kaler and Bessel & Brett has 
been used. 
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distribution. CSP models were based on a scenario where the amplitude of the LF of I 
galaxies, (/)}, increased with redshift as a power law; i.e., (f)*j oc (1 + z) v . The k-corrections 
for the different classes of galaxies were obtained from the spectral energy distributions 
(SEDs) by Coleman, Wu, & Weedman (1980) and Pence (1976); in particular, the SED for 
an average elliptical galaxy was taken to be the SED of the bulge of M31. The e-corrections 
were evaluated from our standard approach to the photometric evolution of galaxies. 

3.1 The Interacting Galaxies 

The I group does not form properly a morphological group of galaxies, it is rather composed 
of pairs of strongly interacting galaxies of different nature. Moreover, the observed physical 
characteristics of the faint galaxies such as: their blue colors, enhanced star formation, 
indicated by [Oil] equivalent widths greater than 20 A (Colless et al. 1994), etc. give 
us information about the nature of the galaxies that intervene in the interactions. These 
characteristics, which are often associated to a system of interacting galaxies, seem to indicate 
that the type of galaxies that are participating in the mergers are gas-rich galaxies (e.g. 
Carlberg & Chariot 1992). This is a point worth noting since we will see below that the 
increase of the characteristic luminosity of the I galaxies, L* T , to a first approximation, is due 
to the starbursts driven by cloud-cloud collision. It will be assumed that the I-galaxy LF 
(LFI) follows a Schechter form. 

The evolution of LFI is due to three factors: (1) a greater number density, (2) a greater 
gas content that results in a stronger starburst, and (3) a smaller average size. The points (2) 
and (3) go in opposite direction, whereas the point (2) increases L* T the point (3) decreases 
it. It is interesting to note that in the absence of starbursts the effects on 5-counts by points 
(1) and (3) cancels out. This might be the case for the i^-counts where most of the light 
comes from a population of old stars. 

Our merger model scenario for the galaxy counts is then the following: suppose that at 
each time a photograph of the Universe is taken and that galaxies are divided into various 
types. Those that show signs of strong interaction are brought to the I group of galaxies. 
As we look back in time we see that the number density (per comoving volume) of I galaxies 
increases, while, by hypothesis, the number density of the other types of galaxies does not 
change. The increase in the total luminosity of I galaxies (oc L* T (j)}) is due to the their greater 
amount of g clS clS well as their increasing number density. 

In this paper the above scenario for the modeling of the B-, i^-counts, and the redshift 
distribution is used and compared with observations. We propose that the LJ of both the 
.B-band and the K-band evolves due to galaxy-galaxy collision and due to the change in 
size, yet a model for the K-counts where the contribution from starbursts is neglected is also 
computed. The galaxy-galaxy collision term, the one which produces an 'active' luminosity 
evolution, contributes with (1 + z) 3 to and it comes from assuming that: (a) L\ is 
proportional to the stellar formation, \1>, driven by the mergers, (b) in a cloud-cloud collision 
\I> oc a 2 , where a is the gas density, and (c) f2 = 1- On the other hand, L} is reduced by 
(1 + z)~ v due to the change in size. Models with different values of rj are calculated. 

A Schechter analytical representation is used for the LFI, with relative freedom in choos- 
ing its a and M* parameters due to the great uncertainty about its shape. Its characteristic 
luminosity, the luminosity at the "knee" of the LF, is assumed to be greater than the corre- 
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sponding value for elliptical and spiral galaxies. The fraction of I galaxies that contributes 
to the LFG is assumed greater than the amount of mergers observed at the present epoch 
(1-2 %) in order to improve our fit to the S-counts. 

3.2 Ingredients 

In Table 1 the values of the parameters of the luminosity functions of the different galaxy 
classes are shown: second, third, and fourth columns. In the fifth, sixth, an seventh columns 
the timescale of the SFR, the B — V and B — K color indices are shown, respectively. These 
later are computed using the photometric evolution code described in §2. The last column 
indicates the kind of SED used for each type of galaxy; in particular, the SED of the bulge 
of M31 is used for E/SO galaxies. 

The assumed values for a, M*, and the mix reflect several facts. First, the values for the 
parameters a and M* of the LF of spiral galaxies are similar to those derived by Loveday 
et al. (1992) for the field LFG. Second, there is still much uncertainty about the value of 
a for early-type galaxies, the rather low value assumed in this paper, as compared with 
those obtained by several authors (Efstathiou et al. 1988; Loveday 1992; Zucca, Pozzetti, & 
Zamorani 1994), fits better the i\~-counts. Third, there is not yet a field LF for early-type 
dwarf galaxies nor for late-type dwarfs, though we use the a and M* values of the LF of 
dE/dSO galaxies and the M* value of the LF of dl galaxies obtained by Sandage, Binggeli, 
& Tammann (1985) for the Virgo cluster; needless to say that there is no reason to expect 
that these assumptions are correct. Fourth, the parameters for the LFI are chosen in such 
a way that: (1) the characteristic luminosity is greater than the one for E/SO galaxies, (2) 
the slope at the faint end equals that of the total luminosity function, and (3) the fraction 
of I galaxies should be small. At the end, our mix takes into account the following facts: 
(1) the uncertainty in the fraction of strongly interacting systems, (2) there are more spiral 
galaxies than elliptical galaxies in the surveys of Efstathiou et al. and Loveday et al., (3) 
there is a significant percentage of dE/dSO galaxies in clusters (Ferguson & Sandage 1991), 
in the poorest cluster, the Leo group, about 46 % are dE/dSO galaxies, and (4) recent 
determinations of the LFG (Eales 1993; Lonsdale & Chokshi 1993) are producing a steeper 
slope at the faint end (~ —1.3). 

TABLE 1 

INPUT PARAMATERS OF THE MODELS 



Type 


a 


M* 


Mix 


(Gyr) 


B- V b 


B — K b 


SED C 


I 


-1.20 


-21.2 


0.05 


oo 


0.49 


3.04 


Sdm 


dl 


-1.40 


-16.2 


0.16 


oo 


0.49 


3.04 


Sdm 


Sa-c 


-0.80 


-20.0 


0.30 


5 


0.70 


3.54 


Scd 


dE/dSO 


-1.35 


-18.0 


0.30 


5.0 


0.70 


3.54 


Sab 


E/SO 


-1.00 


-20.7 


0.19 


0.5 


1.01 


4.17 


M31 



a The star formation rate used to get the color of the galaxies at the present epoch is: tp oc e~*/ T . 

b Two local color indices representatives of the different types of galaxies, computed by our photometric 

evolution of galaxies code. 
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c The spectral energy distributions, SEDs, are from CWW and Pence 1976; in particular, for E/SO galaxies 
we have taken the SED of the bulge of M31. 



The e-corrections are evaluated assuming an exponential decreasing SFR with a varying 
timescale. For dwarf galaxies, spheroidal and irregulars, it appears inappropriate to model 
their light with a solar metallicity population. Furthermore, a choice for the SFR becomes 
difficult; dwarf irregular galaxies are better modeled by a series of star formation bursts 
followed by periods of zero star formation (e.g. Pilyugin 1993). The number and the strength 
of the bursts, and the duration of the quiescent periods are parameters that seem to depend 
on the galaxy. Fortunately, these galaxies make a negligible contribution to the B-, and 
X-counts except for the fainter end (e.g. CSP). The I galaxies have no e-correction by 
hypothesis. 



3.3 Results 

The detection rate and selection effects, according to Yoshii (1993), were evaluated but 
they were not incorporated to our models. The reason is the small effect they have on the 
predicted S-counts. In an isophotal magnitude scheme, Sl = 29 B mag arcsec -2 , with an 
imposed minimum diameter of D min = 2.0", a detection rate of about 75% at B = 26 is 
expected in a (J1 ,A ) = (1,0) scenario. This amounts to a difference of just 0.10 dex at 
B = 26 (see also Yoshii). 

In Figure 3 the number of galaxies versus blue magnitude is plotted for four models. The 
data are from Maddox et al. (1990), Lilly et al. (1991) and Metcalfe et al. (1991). There 
are two non-evolving, NE, models: one with the parameter values from Table 1 (NE) and 
the other one (NE KGB) with the parameter values from Koo, Gronwall, & Bruzual (1993, 
hereafter KGB). The Table 2 from KGB was utilized to derive the Schechter parameters for 
the I (B - V < 0.6), Sa-c (0.6 < B — V < 0.85), and E/SO (B — V > 0.85) classes of 
galaxies, these are: 



(a,M*,<f>*) = < 



(-1.86,-21.6,9.93 x 10~ 5 ) fori 
(—1.13, -21.8, 3.29 x 10- 4 ) for Sa-c 
(-1.38, -20.1, 9.90 x 10- 4 ) for E/SO 



The inclusion of the NE KGB model in Figure 3 is with the purpose of showing that a much 
smaller difference is obtained between a NE model and the observations if the parameters 
of the LFs are moved in the 'correct' direction. The other two models plotted in this Figure 
are merger-driven models for two values of rj. A significant difference is starting to observe 
between each other at bj > 25. Moreover, the maximum difference between the data and 
our rj = 4.0 model amounts just 0.2 dex at bj = 25.0. 

When the LFG is parameterized by a Schechter function, the number density, no, is 
related to the amplitude of the LFG, </>*, by 

n = </>*r(l + a,/3), 

where T is the incomplete gamma function and f3 is the faint luminosity limit in units of 
L*. Our LFG is a composite LF, the sum of five different Schechter LFs, therefore the 
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Figure 3: The number of galaxies versus blue magnitude is plotted for four models: two are 
non-evolving, NE, and the other two are merger-driven, for two values of rj. 
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Figure 4: The logarithm of the LFG is plotted here (dotted line) along with its Schechter fit 
(solid line). The values for the parameters a, 0*, and M* are given in the text. 

parameters a, (ft*, and M* are not well defined. Yet, approximate values for a, <fi*, and M* 
can be obtained by attempting to fit our LFG to a Schechter function. When this is done, 
(see Figure 4) we get the following: a = —1.20, M* = —20.7, and 0* = 1.66 x 10~ 3 Mpc~ 3 . 
It is encouraging to note that the value we get for (/)*, necessary to fit the brighter part of 
the .B-counts, is similar to those derived by several authors (Efstathiou 1988; de Lapparent, 
Geller, & Huchra 1989; Loveday et al. 1992). 

In Figure 5, the number distribution in the X-band is plotted for the same models as 
for the S-band except that here, the NE KGB model is replaced by a merger-driven model 
in which the L* T in the i^-band does not increase due to the startbursts (model denoted 
by II i] = 4.0). The data are from a compilation by Gardner, Cowie, & Wainscoat (1993). 
It is important to point out here the better agreement obtained by those models with a 
rather high value of rj. A redshift distribution is plotted in Figure 6 for two intervals of 
B magnitude. The histograms are redshifts for 125 galaxies taken from Broadhurst et al. 
(1988) in the 20.0 < bj < 21.0 range and for 78 galaxies taken from Colless et al. (1990, 
1993) in the 21.0 < bj < 22.5 interval, this later survey is complete to 4.5 %. Our merger- 
driven models predict a excess of high redshift galaxies, being it greater for models with a 
lower value of r\: the model with r\ = 1.0 predicts that a 27 % of the galaxies have z > 0.7 
in the 21.0 < bj < 22.5 range, clearly inconsistent with the 5 % derived by Colless et al. 
This is interesting since it could be a manner to discriminate between models with number 
luminosity from those with luminosity evolution. The model with rj = 4.0 predicts that 
about 13 % of galaxies have z > 0.7, not too far away from the observed 5 %. 

It appears difficult to reduce the number of high redshift galaxies by changing the param- 
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Figure 5: Number of galaxies against K magnitude for the same models as Figure 3, except 
that the NE KGB model has been replaced by one whose contribution to L* in the if-band 
due to starbursts is neglected (II 77 = 3.0). 

eters of the LF of interacting galaxies, subject to the constraints that: (1) a small percentage 
(< 5%) of galaxies are strongly interacting at the present time, (2) the slope at the faint 
end is that of the total luminosity function, and (3) the characteristic luminosity of this 
population, L*j, is greater than the one of E/SO galaxies. For the sake of the argument, let 
us assume we increase further the value of the parameter 77. One might think this would 
give us a better fit, but the point is that the percentage of galaxies greater than 0.7 does not 
decrease much by increasing the value of 77; for example, for 77 = 10 it is just 12% compared 
to the 13% for 77 = 4!. This point can be explained as follows: as the 77 value goes up the L} 
decreases as 

M; = M;(0) - 2.5(3 - rj)log{\ + z), 

so the number of high-z galaxies decreases (improving the fit) but at the same time the 
number of moderate z galaxies increases due to the (1 + z) v law (worsening the fit). At the 
end, the decreasing of the percentage of galaxies with z > 0.7 is very small. We took the 
value of 77 = 4 as a resonable estimation of the amount of mergers at z = 0.4 and as a not 
too bad fit to the z-distribution. 
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Figure 6: Redshift distributions in two intervals of B magnitude: (a) 20.0 < bj < 21.0, and 
(b) 21.0 < bj < 22.5. Data are from Broadhurst et al. (1988) and Colless et al. (1990; 1993). 
The approximate color relationships B = bj + 0.18 and B = bj + 0.16(bj — r F ) for Broadhurst 
et al. and Colless et al. data, has been used, respectively. The two merger-driven and our 
NE models are plotted. The line code is the same as in Fig. 3. 
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4 CONCLUSIONS 



Differences of up to 0.25 mag in V-K are found when various photometric calibrations are 
used. On the other hand, a smaller difference, 0.1 mag in V-K, is found when we compare the 
results that come from using two reasonable IMFs, the so-called KTG IMF and the Salpeter 
one. This is certainly negligible as far as the B or K number distributions calculations are 
concerned. 

By using the best parameters we have in the literature for the LFs of the different classes 
of galaxies and for the mix of them, we arrive to the conclusion that number and luminosity 
evolution of the galaxy luminosity function are important ingredients in the explanation 
of the excess observed in the S-counts. Our merger-driven models require little evolution 
both in number and luminosity; for example, at z — 0.4 the total number of galaxies has 
just increased by 14 % and the percentage of I galaxies is just 17 % (model with 77 = 4.0). 
Moreover, luminosity evolution is such that a model with 77 = 4.0 predicts that about 13 % of 
galaxies have z > 0.7 in the 21.0 < < 22.5 range. Models with more emphasis in number 
evolution reproduce better the redshift distributions (the 77 = 4.0 model as compared with 
the 77 = 1.0 one), as it should be expected. Despite our relative success of our models, the 
excess of high redshift of galaxies predicted by them indicates us that number and luminosity 
evolution may not be the whole story. 
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Abstract 

A newly developed isochrone synthesis algorithm for the photometric evolution of 
galaxies is described. Two initial mass functions, IMFs, in particular, the recent IMF 
determined by Kroupa, Tout, and Gilmore, three photometric transformations, and a 
1-Gyr-burst star formation rate, SFR, are used to compute the B — V and V — K color 
index evolution. Non-negligible differences are observed among model results. 

In the framework of the galaxy count model by Colin, Schramm, and Peimbert a 
simple merging scenario is considered to account for the excess of galaxies observed in 
the blue band counts. The excess is explained by the number and luminosity evolution 
of a group of galaxies called interacting, I. It is assumed that the number of I galaxies 
increases as (1 + z)' n due to mergers. Moreover, it is proposed that their characteristic 
luminosity increases as (1 + z) 3 due to starbursts driven by galaxy-galaxy collision and 
decreases as (1 + z)~' n due to the change in the size of the galaxies. Not much number 
evolution is needed to account for the excess; for example, a model with r] = 4.0 predicts 
that about 17 % of the galaxies at z = 0.4 are interacting. Number evolution models 
with a rather high value of r] fit better the data; in particular, the model with r] = 4.0 
predicts that about 13 % of the galaxies have z > 0.7 in the 21.0 < ra^ < 22.5 interval, 
this contrasts with the upper bound of 5 % obtained with the sample of 78 galaxies 
by Colless et al. The excess of high redshift galaxies can not be simply explained 
by changing reasonably the parameters of the luminosity function of I galaxies. This 
result could indicate that mergers are not the whole story. Our best-fit model produces 
the following values for the parameters of the local luminosity function of galaxies: 
a = -1.20, M* = -20.7, and <f>* = 1.66 X lO^Mpc" 3 (h = 0.5). 

1 INTRODUCTION 

The excess number density of galaxies observed in the blue band (e.g. Tyson 1988; Metcalfe 
et al. 1991; Lilly, Cowie, & Gardner 1991) over the expected from a non-evolving model has 
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generated great interest. Several possible explanations have been suggested to account for 
it: from those that invoke a nonzero cosmological constant that provides a greater volume 
to accommodate the excess number of galaxies (Yoshii & Takahara 1988; Fukugita et al. 
1990; Yoshii 1993), to those that assume evolution of the luminosity function of galaxies 
(LFG) via an increase in the number density of galaxies due to mergers (Guiderdoni & 
Rocca-Volmerange 1987; Broadhurst, Ellis & Glazebrook 1992; Colin, Schramm, & Peimbert 
1994, hereafter CSP) or via an increase in the characteristic luminosity of the LFG due to 
starbursts driven possibly by galaxy-galaxy collisions (Carlberg & Chariot 1992). Moreover, 
faint redshift surveys by Broadhurst et al. (1988) and Colless et al. (1990, 1993) place 
galaxies quite near, the median redshift is just ~ 0.2 and ~ 0.3, respectively. In addition, 
no galaxies were found with z > 0.7 brighter than bj = 22.5 in the sample by Colless et al. 
This result imposes a strong constraint on evolving models. In these models a non-negligible 
fraction of galaxies with z > 0.7 is expected, indeed larger than the one predicted by a 
non-evolving model. 

With the advent of infrared arrays and infrared CCD detectors it is also possible to 
count galaxies deeply in the A^-band (Cowie et al. 1993; Soifer et al. 1994). Surprisingly, 
a non-evolving model fits relatively well the Accounts. This almost rules out models with 
nonstandard cosmologies, as the ones with a non-zero cosmological constant, since in this 
case they overestimate the number of galaxies. Yet, the observational restrictions admit one 
more explanation, namely a new population of blue galaxies at z > 0.2 that faded or were 
disrupted so as to become invisible at present (Cowie, Songaila, & Hu 1991; Babul & Rees 
1992). 

Local low surface brightness galaxies have similar physical properties to those observed 
in the faint blue excess population (McGaugh 1994). McGaugh suggests that the excess is 
made up of low surface brightness galaxies and his idea rests on a bivariate LFG; that is, a 
LFG that, in addition to type, depends on surface brightness. 

The result that faint blue galaxies (FBGs) are weakly clustered (Efstathiou et al. 1991) 
has been sometimes taken as an evidence against the merger-driven explanation for the 
observed excess of FBGs (see, for instance, Babul & Rees 1992; McGaugh 1994); the low 
value for the angular correlation function at 30" is explained by Efstathiou et al. assuming 
that most FBGs belong to a population that is weakly clustered and intrinsically faint at 
the present epoch. This result is challenged by a recent work by Cole et al. (1994) where 
no evidence is found for the evolution of the comoving correlation length with redshift to 
B = 22, however a more recent paper by Infante & Pritchet (1994) supports the earlier 
results by Efstathiou et al. 

The hypothesis that the observed excess of FBGs is due to an increase in the number 
density of galaxies that are strongly interacting, the merger-driven excess, is retaken here. A 
phenomenological model was developed in CSP to account for three different observational 
restrictions, namely: (1) the i?-counts, (2) the Accounts, and (3) a i?-redshift distribution. 
It was due to the increase in the number density of a type of galaxies, called interacting (I), 
that the observations were fitted. 

This kind of model, like others of similar nature, is supported by a recent work on 
high-resolution imaging of faint blue galaxies (Colless et al. 1994) and by HST studies of 
morphology of high redshift galaxies (Griffiths et al. 1994a, b). Colless et al. find, from a 
sample of 26 galaxies which belong to the redshift survey of Colless et al. (1993), 17 galaxies 
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that have an enhanced star formation rate indicated by [Oil] equivalent widths greater than 
20A. About 30% of these galaxies have companions at projected distances closer than 10 h -1 
kpc. 

This model differs from the one by CSP in the following: (1) a different algorithm for 
the computation of the 'passive' luminosity evolution correction is applied, and (2) although 
the total luminosity of the I galaxies, Lj , here is also not conserved, it is not due to that we 
do not take into account the change in size of the I galaxies, as it occurs in CSP, but due to 
the starsbursts driven by galaxy-galaxy collisions. This is what we call 'active' luminosity 
evolution. The model in this paper also differs from the one by CSP in the assumed values 
for the parameters of the luminosity functions of the different types of galaxies; in particular, 
those values assumed for the characteristic luminosities of Sa-c and E/SO galaxies. Because 
they are fainter in CSP, they produce a better fit for the redshift distribution than the one 
we show in this paper. 

An outline of the paper follows: in §2 we briefly discuss our revised color evolution code 
with an isochrone synthesis algorithm, in §3 our merger-driven model is described, and in §4 
we present our conclusions. 

2 COLOR EVOLUTION OF GALAXIES 

2.1 Model 

Our photometric evolution of galaxies code is fully discussed in CSP therefore here just a brief 
review will be given. Unlike CSP where a standard procedure was used (e.g. Bruzual 1983) 
here an isochrone synthesis algorithm is implemented. The advantages of this algorithm 
have been discussed elsewhere (e.g. Chariot & Bruzual 1991), we just would like to mention 
the usefulness of the algorithm in computing color evolution properties for short (~ 10 8 yr) 
time-scales of the star formation rate. 

2.1.1 Ingredients 

The ingredients necessary in any spectrophotometric evolution of galaxies code are four: (1) 
a library of stellar evolutionary tracks, (2) a library of stellar spectra, (3) an initial mass 
function (IMF), and (4) a star formation rate (SFR). Our compiled tracks are from Schaller 
et al. (1992) for solar metallicity. These tracks are incomplete and therefore we had to 
complete them. Tracks from other authors were used and various useful simplifications were 
made; in particular, the lifetime of the asymptotic giant branch was left as a free parameter. 
This will be remedied as soon as the evolutionary models for the horizontal and asymptotic 
giant branch from the group of Maeder and collaborators are incorporated. As long as we 
are only interested in galaxy color evolution there is no need for having a library of stellar 
spectra, we just need a photometric transformation from the HR theoretical diagram to the 
observational one. The compilation work by Schmidt-Kaler (1982) has been used now to 
pass from one diagram to the other. We still continue to use the work by Johnson (1966) to 
get the rest of the color indexes. 
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The IMF is normalized so that, 



j-m up 

/ m<f)(rn)drn = 1. (1) 

■''miow 

Although a value of 0.1 M & is often used for m; otu , as we do, we should keep in mind that 
the contribution by substellar objects to the integrated infrared light of a galaxy might not 
be negligible. Of course, this depends on the luminosity function of these objects of which 
we know nothing. From chemical evolution considerations Peimbert et al. (1994) limit the 
amount of substellar objects in the solar neighborhood to < 0.02 M & pc~ 3 . The recent IMF 
determination by Kroupa, Tout, and Gilmore (1993), KTG IMF, for the solar vicinity is used, 
as a good guess, for the IMF for a variety of types of galaxies. Nevertheless, the "standard" 
Salpeter IMF, S IMF, with the exponent x = 1.35 is used for comparison purposes. A value 
of 60.0 Mq is taken for m up . 

In view of the great uncertainties about stellar formation histories of galaxies, an expo- 
nential law for the SFR is usually taken to model the color evolution of the different types 
of galaxies. The maximum virtue of this scheme is that almost all present galaxy colors can 
be reproduced by varying the time scale, r, of the SFR: from an elliptical galaxy, r = 0.5, 
to an irregular one, r = oo (e.g. Chariot & Bruzual 1991; CSP). These galaxy colors are, 
of course, metallicity-dependent. This effect is now being incorporated in the latest models 
for spectral evolution of galaxies; in particular, for elliptical galaxies see Worthey (1994 and 
references therein). The evolutionary correction, e-correction, to the counts should not worry 
us too much because as we know from CSP there are only three types of galaxies that play 
an important role on the counts. According to our count-scheme, the luminosity of a galaxy 
belonging to the interacting group of galaxies does not evolve; i.e., it has no e-correction. 
We can be confident that the other two types, E/SO and Sa-c, are being modeled relatively 
well; the first, with a burst of star formation that lasts 1 Gyr, the second, with a constant 
SFR. 



2.1.2 Isochrone synthesis algorithm 

The luminosity of a single stellar population (SSP) is given by 

l AX (t) = 10-°- 4 ( M ^- M ^.°V(™Mm, (2) 

•J Tfl rn i rL 

where m m i n is the lower limit of the IMF and m max is the maximum mass in the isochrone 
with age t. The integrated luminosity for a composite stellar population, with a history of 
star formation given by ip(t), can be obtained with the convolution integral 

L A x(t) = f tP(t-t')l AX (t')dt'. (3) 
Jo 

The isochrone at time t from a complete set of m evolutionary tracks is built as follows: let 
us assume that the tracks are divided into n stages, from zero age main sequence (ZAMS) 
to post-asymptotic giant branch (P-AGB). The indexes i and j number the stages and the 
tracks, respectively, with i = 1, ...,n and j = 1, ...,m. The mass of the star at the stage i is 
given by 

log mi (t) = A hJ log(M j+1 ) + (1 - A hJ ) log(M,), (4) 
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where 

A 8J = lo S**- lo g* . (5) 

logt jti - log t j+lti 

In equation (5) t h i denotes the age of the star of ZAMS-mass Mj at the i th evolutionary 
stage, where 

and 

Mj < m t (t) < M j+1 . 

Mi and Af m denote the lower and upper limit of the IMF, respectively. This procedure 
is performed for each stage (value of i). The HR diagram physical parameters, log L and 
logT e jj, associated to the star of mass ra 4 - are obtained by interpolating the tracks between 
the stars of masses Mj and Mj + \. 

To obtain the integrated properties of this SSP we proceed as Bruzual (1992): the number 
of stars of mass ra 4 -, n(m 8 ), is found by integrating the IMF from m~ to m + , where m~ = 
^Jmi-irrii and m + = ^/m 8 m 8 _|_i. The luminosity, /aa,;, is then computed by assigning Nij = 
Ai } jn(mi) stars to the HR parameters of the star of mass M J+ i at the i th stage, and iV 8J +i = 
(1 — Ai^)n(rrii) to the star of mass Mj at the same i th stage. We finally obtain l&\(t) by 
summing over all i's. 

2.2 Results 

In doing photometric evolution of galaxies special care should be put on the photometric 
transformation. This is stressed in Figure 1 where the B-V and V-K evolution is plotted for 
a 1-Gyr-burst SFR, for three different calibrations. The dotted line uses a compilation work 
by Schmidt-Kaler (1982) and Johnson (1966), the solid line utilizes one from Schmidt-Kaler 
(1982) and Bessel and Brett (1988), and the dashed line comes from a compiled calibration 
work by Bruzual and Chariot (1993). The greatest difference among them amounts to ~ 0.25 
mag in V-K at around 1.7 X 10 7 yr; this difference is due to the V-K values adopted for the 
giant luminosity class at the low-temperature regime. The peak at 1.7 X 10 7 yr is produced by 
the supergiant red branch contribution. Hereafter the calibration by Schmidt-Kaler (1982) 
and Bessel and Brett (1988) will be used. In Figure 2 we have plotted the B-V and V-K color 
evolution for a 1-Gyr-burst SFR, in this case for two different IMFs: the KTG IMF (solid 
line) and the Salpeter one (dotted line). The V-K color is slightly bluer when computed 
with the KTG IMF, the greatest difference at ~ 10 8 j/r amounts to 0.1 mag. A much smaller 
difference is observed for the case of B-V. 



3 A Merger Model for Galaxy Counts 

This paragraph describes briefly the fi = 1 merger model for galaxy counts of CSP. Galaxies 
were divided into five different types: E/SO, dE/dSO (spheroidal dwarfs), Sa-c (spirals), dl 
(dwarf irregulars), and I (interacting galaxies). A Schechter analytical representation for the 
galaxy luminosity function was used for all types. It was shown that a gaussian representation 
for the luminosity function (LF) of the brighter galaxies provided a poor fit to the redshift 
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Figure 1: Evolution of the color indices B-V and V-K for a lGyr-burst SFR with three 
different photometric calibrations: Bruzal & Chariot (dashed line), Schmidth-Kaler and 
Bessel & Brett (solid line), and Schmidth-Kaler and Johnson (dotted line). 
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Figure 2: As in Fig. 1 we have plotted the evolution of the color indices B-V and V-K for 
a 1-Gyr-burst SFR. The two curves come from two different IMFs: KTG IMF (solid line) 
and Salpeter IMF (dotted line). The calibration from Smichdt-Kaler and Bessel & Brett has 
been used. 
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distribution. CSP models were based on a scenario where the amplitude of the LF of I 
galaxies, <f>* 1} increased with redshift as a power law; i.e., <f>*j oc (1 + z)' n . The k-corrections 
for the different classes of galaxies were obtained from the spectral energy distributions 
(SEDs) by Coleman, Wu, & Weedman (f980) and Pence (1976) ; in particular, the SED for 
an average elliptical galaxy was taken to be the SED of the bulge of M31. The e-corrections 
were evaluated from our standard approach to the photometric evolution of galaxies. 

3.1 The Interacting Galaxies 

The I group does not form properly a morphological group of galaxies, it is rather composed 
of pairs of strongly interacting galaxies of different nature. Moreover, the observed physical 
characteristics of the faint galaxies such as: their blue colors, enhanced star formation, 
indicated by [Oil] equivalent widths greater than 20 A (Colless et al. 1994), etc. give 
us information about the nature of the galaxies that intervene in the interactions. These 
characteristics, which are often associated to a system of interacting galaxies, seem to indicate 
that the type of galaxies that are participating in the mergers are gas-rich galaxies (e.g. 
Carlberg & Chariot 1992). This is a point worth noting since we will see below that the 
increase of the characteristic luminosity of the I galaxies, L* I} to a first approximation, is due 
to the starbursts driven by cloud-cloud collision. It will be assumed that the I-galaxy LF 
(LFI) follows a Schechter form. 

The evolution of LFI is due to three factors: (1) a greater number density, (2) a greater 
gas content that results in a stronger starburst, and (3) a smaller average size. The points (2) 
and (3) go in opposite direction, whereas the point (2) increases L\ the point (3) decreases 
it. It is interesting to note that in the absence of starbursts the effects on i?-counts by points 
(1) and (3) cancels out. This might be the case for the Accounts where most of the light 
comes from a population of old stars. 

Our merger model scenario for the galaxy counts is then the following: suppose that at 
each time a photograph of the Universe is taken and that galaxies are divided into various 
types. Those that show signs of strong interaction are brought to the I group of galaxies. 
As we look back in time we see that the number density (per comoving volume) of I galaxies 
increases, while, by hypothesis, the number density of the other types of galaxies does not 
change. The increase in the total luminosity of I galaxies (oc L}^}) is due to the their greater 
amount of e well as their increasing number density. 

In this paper the above scenario for the modeling of the B- } Accounts, and the redshift 
distribution is used and compared with observations. We propose that the L} of both the 
i?-band and the A^-band evolves due to galaxy-galaxy collision and due to the change in 
size, yet a model for the Accounts where the contribution from starbursts is neglected is also 
computed. The galaxy-galaxy collision term, the one which produces an 'active' luminosity 
evolution, contributes with (1 + z) 3 to L\ and it comes from assuming that: (a) L\ is 
proportional to the stellar formation, ^ , driven by the mergers, (b) in a cloud-cloud collision 
^ oc <t 2 , where a is the gas density, and (c) fi = 1. On the other hand, L} is reduced by 
(1 + z)~' n due to the change in size. Models with different values of rj are calculated. 

A Schechter analytical representation is used for the LFI, with relative freedom in choos- 
ing its a and M* parameters due to the great uncertainty about its shape. Its characteristic 
luminosity, the luminosity at the "knee" of the LF, is assumed to be greater than the corre- 
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sponding value for elliptical and spiral galaxies. The fraction of I galaxies that contributes 
to the LFG is assumed greater than the amount of mergers observed at the present epoch 
(1-2 %) in order to improve our fit to the i?-counts. 

3.2 Ingredients 

In Table 1 the values of the parameters of the luminosity functions of the different galaxy 
classes are shown: second, third, and fourth columns. In the fifth, sixth, an seventh columns 
the timescale of the SFR, the B — V and B — K color indices are shown, respectively. These 
later are computed using the photometric evolution code described in §2. The last column 
indicates the kind of SED used for each type of galaxy; in particular, the SED of the bulge 
of M31 is used for E/SO galaxies. 

The assumed values for a, Af*, and the mix reflect several facts. First, the values for the 
parameters a and M* of the LF of spiral galaxies are similar to those derived by Loveday 
et al. (1992) for the field LFG. Second, there is still much uncertainty about the value of 
a for early-type galaxies, the rather low value assumed in this paper, as compared with 
those obtained by several authors (Efstathiou et al. 1988; Loveday 1992; Zucca, Pozzetti, & 
Zamorani 1994), fits better the Accounts. Third, there is not yet a field LF for early-type 
dwarf galaxies nor for late-type dwarfs, though we use the a and M* values of the LF of 
dE/dSO galaxies and the M* value of the LF of dl galaxies obtained by Sandage, Binggeli, 
& Tammann (1985) for the Virgo cluster; needless to say that there is no reason to expect 
that these assumptions are correct. Fourth, the parameters for the LFI are chosen in such 
a way that: (1) the characteristic luminosity is greater than the one for E/SO galaxies, (2) 
the slope at the faint end equals that of the total luminosity function, and (3) the fraction 
of I galaxies should be small. At the end, our mix takes into account the following facts: 
(1) the uncertainty in the fraction of strongly interacting systems, (2) there are more spiral 
galaxies than elliptical galaxies in the surveys of Efstathiou et al. and Loveday et al., (3) 
there is a significant percentage of dE/dSO galaxies in clusters (Ferguson & Sandage 1991), 
in the poorest cluster, the Leo group, about 46 % are dE/dSO galaxies, and (4) recent 
determinations of the LFG (Eales 1993; Lonsdale & Chokshi 1993) are producing a steeper 
slope at the faint end (~ —1.3). 

TABLE 1 

INPUT PARAMATERS OF THE MODELS 



Type 


a 


M* 


Mix 


(Gyr) 


B - V b 


B - K b 


SED C 


I 


-1.20 


-21.2 


0.05 


oo 


0.49 


3.04 


Sdm 


dl 


-1.40 


-16.2 


0.16 


oo 


0.49 


3.04 


Sdm 


Sa-c 


-0.80 


-20.0 


0.30 


5 


0.70 


3.54 


Scd 


dE/dSO 


-1.35 


-18.0 


0.30 


5.0 


0.70 


3.54 


Sab 


E/SO 


-1.00 


-20.7 


0.19 


0.5 


1.01 


4.17 


M31 



a The star formation rate used to get the color of the galaxies at the present epoch is: ip oc e _ */ r . 

b Two local color indices representatives of the different types of galaxies, computed by our photometric 

evolution of galaxies code. 
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c The spectral energy distributions, SEDs, are from CWW and Pence 1976; in particular, for E/SO galaxies 
we have taken the SED of the bulge of M31. 

The e-corrections are evaluated assuming an exponential decreasing SFR with a varying 
timescale. For dwarf galaxies, spheroidal and irregulars, it appears inappropriate to model 
their light with a solar metallicity population. Furthermore, a choice for the SFR becomes 
difficult; dwarf irregular galaxies are better modeled by a series of star formation bursts 
followed by periods of zero star formation (e.g. Pilyugin 1993). The number and the strength 
of the bursts, and the duration of the quiescent periods are parameters that seem to depend 
on the galaxy. Fortunately, these galaxies make a negligible contribution to the and 
Accounts except for the fainter end (e.g. CSP). The I galaxies have no e-correction by 
hypothesis. 

3.3 Results 

The detection rate and selection effects, according to Yoshii (1993), were evaluated but 
they were not incorporated to our models. The reason is the small effect they have on the 
predicted i?-counts. In an isophotal magnitude scheme, Sl = 29 B mag arcsec -2 , with an 
imposed minimum diameter of D m i n = 2.0", a detection rate of about 75% at B = 26 is 
expected in a (fi ,Ao) = (1,0) scenario. This amounts to a difference of just 0.10 dex at 
B = 26 (see also Yoshii). 

In Figure 3 the number of galaxies versus blue magnitude is plotted for four models. The 
data are from Maddox et al. (1990), Lilly et al. (1991) and Metcalfe et al. (1991). There 
are two non-evolving, NE, models: one with the parameter values from Table 1 (NE) and 
the other one (NE KGB) with the parameter values from Koo, Gronwall, & Bruzual (1993, 
hereafter KGB). The Table 2 from KGB was utilized to derive the Schechter parameters for 
the I (B - V < 0.6), Sa-c (0.6 < B - V < 0.85), and E/SO (B - V > 0.85) classes of 
galaxies, these are: 

C (-1.86,-21.6,9.93 x 10" 5 ) fori 
(a, M* ,</>*)= I ( — 1.13, —21.8,3.29 x 10" 4 ) for Sa-c 
( (-1.38,-20.1,9.90 x 10" 4 ) for E/SO 

The inclusion of the NE KGB model in Figure 3 is with the purpose of showing that a much 
smaller difference is obtained between a NE model and the observations if the parameters 
of the LFs are moved in the 'correct' direction. The other two models plotted in this Figure 
are merger-driven models for two values of rj. A significant difference is starting to observe 
between each other at bj > 25. Moreover, the maximum difference between the data and 
our j] = 4.0 model amounts just 0.2 dex at bj = 25.0. 

When the LFG is parameterized by a Schechter function, the number density, n , is 
related to the amplitude of the LFG, </>*, by 

n = <f)*T(l + a, /3), 

where T is the incomplete gamma function and f3 is the faint luminosity limit in units of 
L*. Our LFG is a composite LF, the sum of five different Schechter LFs, therefore the 
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Figure 3: The number ol galaxies versus blue magnitude is plotted for four models: two are 
non-evolving, NE, and the other two are merger- driven, for two values of r\. 
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Figure 4: The logarithm of the LFG is plotted here (dotted line) along with its Schechter fit 
(solid line). The values for the parameters a, </>*, and M* are given in the text. 

parameters a, </>*, and M* are not well defined. Yet, approximate values for a, </>*, and Af* 
can be obtained by attempting to fit our LFG to a Schechter function. When this is done, 
(see Figure 4) we get the following: a = -1 .20, M* = -20.7, and <j>* = 1 .66 X 10~ 3 Mpc~ 3 . 
It is encouraging to note that the value we get for </>*, necessary to fit the brighter part of 
the i?-counts, is similar to those derived by several authors (Efstathiou f988; de Lapparent, 
Geller, & Huchra f989; Loveday et al. 1992). 

In Figure 5, the number distribution in the A^-band is plotted for the same models as 
for the i?-band except that here, the NE KGB model is replaced by a merger-driven model 
in which the L\ in the A^-band does not increase due to the startbursts (model denoted 
by II j] = 4.0). The data are from a compilation by Gardner, Cowie, & Wainscoat (1993). 
It is important to point out here the better agreement obtained by those models with a 
rather high value of rj. A redshift distribution is plotted in Figure 6 for two intervals of 
B magnitude. The histograms are redshifts for 125 galaxies taken from Broadhurst et al. 
(1988) in the 20.0 < bj < 21.0 range and for 78 galaxies taken from Colless et al. (1990, 
1993) in the 21.0 < bj < 22.5 interval, this later survey is complete to 4.5 %. Our merger- 
driven models predict a excess of high redshift galaxies, being it greater for models with a 
lower value of rj: the model with rj = f.O predicts that a 27 % of the galaxies have z > 0.7 
in the 21.0 < bj < 22.5 range, clearly inconsistent with the 5 % derived by Colless et al. 
This is interesting since it could be a manner to discriminate between models with number 
luminosity from those with luminosity evolution. The model with rj = 4.0 predicts that 
about f3 % of galaxies have z > 0.7, not too far away from the observed 5 %. 

It appears difficult to reduce the number of high redshift galaxies by changing the param- 
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Figure 5: Number of galaxies against K magnitude for the same models as Figure 3, except 
that the NE KGB model has been replaced by one whose contribution to L* in the A^-band 
due to starbursts is neglected (II rj = 3.0). 

eters of the LF of interacting galaxies, subject to the constraints that: (1) a small percentage 
(< 5%) of galaxies are strongly interacting at the present time, (2) the slope at the faint 
end is that of the total luminosity function, and (3) the characteristic luminosity of this 
population, L* I} is greater than the one of E/SO galaxies. For the sake of the argument, let 
us assume we increase further the value of the parameter rj. One might think this would 
give us a better fit, but the point is that the percentage of galaxies greater than 0.7 does not 
decrease much by increasing the value of rj; for example, for 77 = 10 it is just 12% compared 
to the 13% for rj = 4!. This point can be explained as follows: as the rj value goes up the L} 
decreases as 

m; = m;(o) - 2.5(3 - rj)iog(i + z), 

so the number of high-z galaxies decreases (improving the fit) but at the same time the 
number of moderate z galaxies increases due to the (1 + z)' n law (worsening the fit). At the 
end, the decreasing of the percentage of galaxies with z > 0.7 is very small. We took the 
value of j] = 4 as a resonable estimation of the amount of mergers at z = 0.4 and as a not 
too bad fit to the z-distribution. 
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Figure 6: Redshift distributions in two intervals of B magnitude: (a) 20.0 < bj < 21.0, and 
(b) 21.0 < bj < 22.5. Data are from Broadhurst et al. (1988) and Colless et al. (1990; 1993). 
The approximate color relationships B = &/ + 0.18 and B = bj + 0.16(6/ — rp) for Broadhurst 
et al. and Colless et al. data, has been used, respectively. The two merger-driven and our 
NE models are plotted. The line code is the same as in Fig. 3. 
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4 CONCLUSIONS 



Differences of up to 0.25 mag in V-K are found when various photometric calibrations are 
used. On the other hand, a smaller difference, 0.1 mag in V-K, is found when we compare the 
results that come from using two reasonable IMFs, the so-called KTG IMF and the Salpeter 
one. This is certainly negligible as far as the B or K number distributions calculations are 
concerned. 

By using the best parameters we have in the literature for the LFs of the different classes 
of galaxies and for the mix of them, we arrive to the conclusion that number and luminosity 
evolution of the galaxy luminosity function are important ingredients in the explanation 
of the excess observed in the i?-counts. Our merger-driven models require little evolution 
both in number and luminosity; for example, at z = 0.4 the total number of galaxies has 
just increased by 14 % and the percentage of I galaxies is just 17 % (model with rj = 4.0). 
Moreover, luminosity evolution is such that a model with rj = 4.0 predicts that about 13 % of 
galaxies have z > 0.7 in the 21.0 < < 22.5 range. Models with more emphasis in number 
evolution reproduce better the redshift distributions (the rj = 4.0 model as compared with 
the j] = 1.0 one), as it should be expected. Despite our relative success of our models, the 
excess of high redshift of galaxies predicted by them indicates us that number and luminosity 
evolution may not be the whole story. 
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